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Abstract. The dynamics of a class of zero-range processes exhibiting a condensation 
transition in the stationary state is studied. The system evolves in time starting from a 
random disordered initial condition. The analytical study of the large-time behaviour 
of the system in its mean-field geometry provides a guide for the numerical study of 
the one-dimensional version of the model. Most qualitative features of the mean-field 
case are still present in the one-dimensional system, both in the condensed phase and 
at criticality. In particular the scaling analysis, valid for the mean-field system at 
large time and for large values of the site occupancy, still holds in one dimension. The 
dynamical exponent z, characteristic of the growth of the condensate, is changed from 
its mean-field value 2 to 3. In presence of a bias, the mean-field value z = 2 is recovered. 
The dynamical exponent z c , characteristic of the growth of critical fluctuations, is 
changed from its mean-field value 2 to a larger value, z c ~ 5. In presence of a bias, 
z c — 3. 
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1. Introduction 

In a recent work, Kafri et al. [1] gave a criterion for the existence of phase separation in 
driven one-dimensional systems. Discussions on this theme have been the subject of a 
long series of investigations [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]. The criterion 
relies on a mapping, in the stationary state, of the system onto a zero-range process. 

For example, consider a two-species model, where each site of a one- dimensional 
ring is either occupied by a (+) particle, or by a (— ) particle, or is vacant (0). Positive 
particles are driven to the right, negative particles to the left, and they exchange 
positions when meeting. Define domains as uninterrupted sequences of (+) or (— ) 
particles. Then identify the zeros with boxes, the particles with balls, and a domain 
with the content of the box attached to the zero to its left, say. Then the system of (±) 
particles can be viewed as a system of balls in boxes, where the rate at which a ball 
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leaves a box is given by the current out of a domain. This process is known as a zero- 
range process (ZRP) [17, 18], more precisely defined below. Note that such a mapping 
of the particle system to the corresponding ZRP gives a coarse-grained description of 
the former, because the spatial structure inside domains is lost. Once the hopping rate 
dk out of a box containing k balls is known, the analysis of the stationary state of the 
ZRP is straightforward. Depending on the behaviour of dj, at large k, it is easy to show 
that there is, or there is not, a condensation transition, where a macroscopic fraction of 
balls is contained in one box, the condensate. This corresponds in the original particle 
system to a condensed phase where a single domain contains a macroscopic fraction of 
the particles. 

This is, in essence, the analysis given in [1]. Examples of particle systems with 
two species satisfying the criterion for the existence of phase separation have been 
subsequently given in [19]. 

These considerations concern the stationary state. We are thus naturally lead 
to address the question of the dynamics, both for the underlying original microscopic 
particle models, and for their coarse-grained counterpart, the ZRP. 

For the latter, the agenda is clear. The question is to determine the temporal 
evolution of the system starting from a random initial condition, once given the hopping 
rate out of a box. The present paper is devoted to the study of the dynamics of a class of 
ZRP giving rise to a condensation transition in their stationary state. We will proceed 
in two steps. We study first the system in its mean-field geometry, where all boxes are 
connected. This gives then a framework for the study of the one-dimensional system. 

The dynamics of the microscopic models defined in [19], or more generally of two- 
species models satisfying the criterion for phase separation, is more subtle because it 
involves two (a priori different) scales of time: the intra-domain dynamics, corresponding 
to the equilibration inside a domain, and the inter-domain dynamics, corresponding to 
the coarsening of domains. One can therefore wonder whether during its temporal 
evolution the particle system is equivalent to a ZRP with constant hopping rate. This 
study is postponed to another publication [20]. 

The aim of this paper is therefore two-fold. Firstly it addresses the question of 
the dynamics of zero-range processes per se. Indeed, ZRP are sufficiently ubiquitous in 
the field of stochastic processes to warrant the interest of such a study. For instance, 
the analysis of the zeta urn model [21, 22] is surprisingly close to that presented in the 
next sections. Secondly it provides a template for the study of the dynamics of particle 
systems, whenever these can be mapped, at least at a coarse-grained level, to a ZRP, as 
explained above. 
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2.1. Mean- field geometry and equilibrium: general formalism 

The zero-range process studied hereafter is defined as follows. Consider a one- 
dimensional lattice of M sites (or boxes), % = 1, . . . , M . In each box we have JVj 
indistinguishable particles such that 

M 

J2 N i = N - 

i=l 

The dynamics of the system is given by the rates at which a particle leaves the departure 
box d, containing Nd = k particles, and is put in the arrival box a containing N a = I 
particles. The hopping rate, denoted by W k i, a priori depends both on k and I. 

A configuration of the system is specified by the occupation numbers Ni(t), hence 
a complete knowledge of its dynamics involves the determination of 7- > (A 1) A 2 , . . . , Nm), 
the probability of finding the system in a given configuration at time t. Hereafter we 
will only focus our attention on the probability of finding k particles in the generic box 
i = l, 

fk(t) = V(N 1 (t) = k), 

that is, the marginal distribution of the occupation number Ni(t) of this box. 
Conservation of probability and of density yields 

oo 

E/*(0 =1 » (!) 

fc=0 

oo 

E */*(*)= ft (2) 

k=l 

where we have taken the thermodynamic limit iV — > oo, M — > oo, with fixed density 
p = N/M. 

Suppose now that the system has the mean- field geometry, i.e., all the boxes are 
connected. In this case the temporal evolution of the occupation probability fk(t) is 
explicitly given by the master equation 
df k (t) 



dt 
d/o(t) 
dt 



— Pfc+l fk+l + Afc_i fk-1 — {fJ>k + (k > 1), 

— A*i fi — Aq/o, 



where 



^=£W fc ,/,, (3) 

A fc =EWi fc /,, (4) 

are respectively the rates at which a particle leaves box number 1 (A^ = k — * iVi = fc— 1), 
or arrives in this box (A^ = fc — > Ni = k + 1). This is the master equation for a random 
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walk for N 1: i.e., on the positive integers k = 0,1, The equation for f is special 

because one cannot select an empty box as a departure box, nor can Ni be negative. 
The conservation law dJ2k fk/dt = 0, valid at all times, implies the sum rule 

oo oo 

Yl f^k fk = Yl ^ fk- (5) 

fe=l fc=0 

In the stationary state (fk = 0), we find 

/fc+l,eq _ Afc 
/fc,eq A*fe+1 

which expresses the detailed balance condition at equilibrium, yielding 
_ Ap . . . Afc-i 

Jfc,eq J0,cqj 
/i! . . . /i fc 

where /o, eq is fixed by normalisation. 

Hereafter we are interested in the case where the rate Wki only depends on the 
occupation of the departure box. We denote the hopping rate out of the departure 
box by Wki = dk- (Alternatively one can consider cases where Wki only depends on the 
occupation of the arrival box, as for the zeta urn model [21, 22]: Wki — o-i is the hopping 
rate into the arrival box; see below.) We thus have 

oo 

Hk = 4, A fc = J2 d ifi = dt- (6) 
i=i 

Note that (5) is automatically satisfied. We obtain 

_ Pk z k 

Jk,eq v00 k , {I) 

where z = J2 dji m = d cq , and 
Po = 1- 

It is important to emphasize here that the results above, equations (7) and (8), 
also hold for the equivalent one-dimensional ZRP in the stationary state. Indeed 
a well-known property of ZRP states that in the stationary state the probability 
V(NiN 2 , . . . , N M ) is given by a product measure [17, 18]. Explicitly V(N 1 N 2 , . . . , N M ), 
with Ni = N, factorizes into the product /(iVi)/(A^) • • • /(-/Vm), where each factor 
f(Ni = k) is equal to fk, eq as given by (7) and (8). For an infinite system this 
factorization property holds at all time in the mean field geometry, by essence. 

2. 2. Class of ZRP with a condensation transition 

Let us apply this formalism to the case where, at large k 

d k w 1 + ^ + • • • ■ (9) 
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Then it is easy to see that, in the same regime, p k ~ k~ b . Hence the series Y,kLoPk zk 
has a radius of convergence z c — 1. The unknown parameter z = d cq is determined by 
using (2), which yields 

£*/M, = §S?^pr = P- do) 

fc=l 2^k=oPkZ 

This expression is monotonously increasing with z. When z reaches the value z c = 1, 
the right side is converging only if b > 2. Therefore we have the following two cases. 

• If b < 2, equation (10) possesses a solution z(p) for any value of p. The occupation 
probabilities are exponentially decaying as 

hn= ^ kZ l 7k - k ~ b ^(-j) ( fluid )' ( n ) 

where the correlation length £ = |lnz| _1 . The system is in a fluid phase. 

• If b > 2, there exists a critical value of the density, where the correlation length 
diverges, attained when z = z c = 1: 

* = W?- < 12 > 

This is the equation of a transition line. In other terms, if p < p c , then (10) has a solution 
z(p), in which case the system is in the fluid phase described above (see equation (11)). 
Otherwise if p > p c , then necessarily z = d cq = 1. The system is composed of a critical 
part and a condensate. The critical part is described by the algebraically decaying 
distribution 

f k , eq = — ^ k~ b (critical) (13) 

contributing to p c , as shown by (12). The condensate, contributing to the remaining 
p — p c , corresponds to particles sitting in a single box. At criticality (p = p c ) the 
condensate disappears and (13) still holds. 

Let us illustrate the above formalism for the explicit choice of hopping rate 

d k = 1 + p (14) 
which is used in the numerical studies of the forthcoming sections. We have 

/ fe e q = (&-l) r( / ? r( f +1 ? ~ (6 - l)r(6) k~\ (15) 

J , q v ; T{k + b+ 1) fc^oo v ; v ; v ; 

and 

Pc = ^ (16) 

which is the equation of the critical line in the b — p plane. It is simple to check using 
(15) that we have consistently d cq = 1, both at criticality and in the condensed phase. 
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From now on, we consider a hopping rate d k out of the departure box of the form (9), i.e., 
such as dk ~ 1 + b/k at large k, and for the sake of analytical or numerical illustrations 
we will choose the explicit form (14), unless specified. 

For an infinite system in the mean-field geometry, the knowledge of the occupation 
probabilities fk(t) provides a complete description of its dynamics. We consider a system 
with Poissonian initial distribution of occupation probabilities, 

A(0) = e"^, 

i.e., such that initially particles are distributed at random amongst boxes. The temporal 
evolution of the occupation probabilities is given by the master equation 
df k (t) 



dt 
dfo(t) 



= d k+1 f k+1 + dtfk-i ~ (dk + d t )fk (k>l), (17) 
— di fi — dtfo- 



dt 

This set of equations is non linear because d t = J^difi is itself a function of the f k (t). 
Hence there is no explicit solution of these equations in closed form. Yet one can extract 
from them an analytical description of the dynamics of the system at long times, both 
in the condensed phase, and at criticality. The structure of the reasoning borrows to 
former studies on urn models [21, 22]. (For a review, see [23].) 

As we show below, there exists two different regimes in the evolution of the system, 
both in the condensed phase or at criticality, which we study successively. 

3.1. Condensed phase (p > p c ) 
Since d eq = 1, we set, for large times, 

dt^l + Aet, (18) 
where the small time scale e t is to be determined. As shown below, e t ~ t~*. 



Regime I: k fixed, t large. For t large enough, boxes empty (dk) faster than they fill 
(dt). In this regime there is convergence to equilibrium, hence we set 

fk(t) « /fc,eq(l +v k e t ), (19) 

with / fc eq given by (13), and where the v k are unknown. This expression carried into 
(17) yields the stationary equation f k = 0, because the derivative proportional to 
e t , is negligible compared to the right-hand side. We thus obtain the detailed balance 
condition: 

/fc+i,cq 1 + Vk+i £t _ 1 + As t 

fk,eq 1 + Vk£t dk+1 
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Using (13) and (8), we obtain, at leading order in e t , Vk+i — v k = A, and finally 

v k = v + kA. (20) 
At this stage, t>o and the amplitude A are still to be determined. 

Regime II: k and t are simultaneously large. This is the scaling regime, with scaling 
variable u — ket- We look for a similarity solution of (17) of the form 

f k (t)*e 2 t g(u). (21) 

Following the treatment of [21] we obtain for g(u) the linear differential equation 

/(«) + {\ - A + -\ g'(u) + (l - A) g(u) = 0. 



f 
Jo 



This is precisely the differential equation found in [21, 22]. Moreover this method leads 
to e t ~ [21]. The amplitude A can be determined by the fact that the equation 
has an acceptable solution g(u) vanishing as u — > and u — > oo [21]. It is a universal 
quantity, only depending on the value of b. The normalisation of the solution is fixed 
by the sum rule (1) yielding [21] 

duug{u) = p — p c . (22) 

The sum rule (2) leads to [21] 

roo 

v + Ap c = - dug(u). (23) 
Jo 

The differential equation above has no closed form solution. However further information 
on the form of the solution g{u) can be found in [21, 22]. In particular g{u)/(p — p c ) 
is universal in the sense that it only depends on the parameter b (see the discussion 
below) . 

An intuitive description of the dynamics of the condensate in the scaling regime is 
as follows. The typical occupancy k cond of the boxes making the condensate scales as 
ts. The total number of particles in the condensate is equal to M(p — p c ), the remaining 
M p c lying in the fluid. Therefore the number of boxes belonging to the condensate 
scales as M(p — p c )t~\. 

3.2. Criticality (p = p c ) 

The analysis follows closely that done in [22]. We just sketch the method here. 
We set 

dt w 1 + Ae t , 

with e t = t _UJ , where the exponent u> is to be determined, and we consider the same 
two regimes as above. In regime I, we still set (19) for fk(t). The reasoning leading to 
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the relationship v k = v + k A (see (20)) is still valid here. In regime II, we look for a 
similarity solution to (17) of the form 

fk(t) ~ fk,eq9c(u) u = ht~*. (24) 

The sum rules (1) and (2) lead respectively to the following equations, provided that 
b > 3, 

v + A Pc = 0, (25) 

t~ u (v oPc + Ap c ) = r^' 2 {b - l)r(6) / &uv}- b {l - g c (u)), (26) 

Jo 

where, using (15), /i c = J^k 2 fk, cq = (b + l)(b — 2) _1 (6 — 3) _1 . Equation (26) fixes the 
value of uj: 

u = {b-2)/2. (27) 
The differential equation obeyed by g c (u) is obtained by carrying (24) into (17). It reads 

sTM + - b ~) 9' c {u) = o, 

the solution of which is, with g c (0) = 1, 

o— b pc<j 

3cW^ dyy b e~y 2 /\ (28) 

1 { — ) Ju 

The fall-off of g c (u) for u ^> 1 is very fast: g c (u) ~ exp(— m 2 /4), hence fk(t) ~ 
exp(— k 2 /4t). We finally obtain 

, (6-l)r(6) r°° , , hl . (6-2)(6-3)„ /6\ 
/j, c - pi Jo 6-1 V 2 / 

5.5. Universality 

To summarize, for any hopping rate of the form c4 1 + b/k, the scaling functions, 
(7(m) in the condensed phase (more precisely: g(u)/(p — p c )), and g c {u) at criticality, 
are universal. In both cases the scaling variable is u — kt~ x l 2 . The critical density p c , 
and, as a consequence, any quantity depending on p c , such as the amplitude Vq, are non 
universal, with values depending on the precise definition of d^. As noted above, the 
amplitude A is a universal quantity in the condensed phase. 

Universality can be checked by comparing the results obtained by a numerical 
integration of the master equation (17), for different choices of the rate d^. Besides the 
form d k — 1 + b/k of eq. (14) we investigated the following two other expressions 

d k = l + ± + ± (29) 



1 



b 



d k = ( 1 + ^ ) • (30) 
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Figure 1. Mean-field geometry: normalised scaling function ug(u)/(p — p c ) 
in the condensed phase, for times equal to 10 3 , 10 4 , and 10 5 . (b = 4, 
p = 20 Pc = 10.) 

Figure 1 depicts the normalized scaling function ug(u)/(p — p c ) in the condensed 
phase, for times equal to 10 3 , 10 4 , and 10 5 . It was obtained with the choice of rate (14), 
with 6 = 4, and for a density p = 20p c . The curves obtained with the choice (29) and 
c = 10, or with (30), are hardly distinguishable from the former one. Figure 2 depicts 
g c {u) at criticality, for the same range of times, with the choice (14), and 6 = 4. Again, 
the curves obtained with the choice (29), with c = 10, or with (30), are indistinguishable 
from the former one. The limiting curve of figure 2 coincides exactly with the form (28). 

One also finds, for p > p c and b = 4, A ~ 1.9, and / dug(u)/(p — p c ) ~ 0.336, 
which are universal quantities. Correspondingly, for dk given by (14), and p = 20p c , 
v ~ —4.15. 

It is remarkable that the universality of the scaling functions g(u) and g c {u) extends 
beyond the case of the ZRP under study. Indeed, for the zeta urn model, defined 
by a rate Wki only depending on the occupation I of the arrival box, much of the 
formalism above is the same [21, 22]. An intuitive explanation of this fact comes from 
the comparison of the rates of the two models, at large values of time and occupation 
numbers. Consider the rate pk and defined in equations (3) and (4). For the ZRP 
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Figure 2. Mean-field geometry: scaling function g c (u) at criticality for times 
equal to 10 3 , 10 4 , and 10 5 . (b = 4, p = p c = 0.5.) 

studied here these rates are given by (cf. (6)) 

b 

fi k = d k ^l + -, X k = Y,difi = dt-l + Ae t , (ZRP) 

For the zeta urn model (see [21, 22] for detailed definitions), we have 

fi k = 1, A fc ^1 - j^j (1 + Ae t ) (zeta), 

where (3 is the inverse temperature, playing the role of the parameter b in the model. 
Simple inspection of the last two equations demonstrates the formal analogy between 
the two models. This analogy is still stronger if one chooses the rate d k given by (30). 
Then, at equilibrium, the formalisms for the corresponding ZRP and for the zeta urn 
model coincide. 

4. ZRP: dynamics in one dimension 

The framework of analytic results obtained in the mean-field case provides a guide 
for the numerical study of the one-dimensional system. It turns out that most of the 
features of the former survive in one dimension. In particular the two regimes of time 
described in the previous section are still present in the one-dimensional case. The main 
difference lies in the values of the dynamical exponent z characteristic of the growth of 
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the condensate, and of the critical dynamical exponent z c characteristic of the growth of 
critical fluctuations. In the condensed phase (p > p c ), numerical simulations, presented 
hereafter, give evidence for z — 3, instead of z = 2 in the mean-field geometry. At 
criticality z c ~ 5 instead of z c = 2, in the mean-field geometry. 

We consider a system on a ring. Sites have multiple occupancy, given by the 
occupation numbers iVj(i), where % = 1, . . . ,M. At each time step a site % is chosen 
at random and, if not empty, a particle of this site is transferred to one of the two 
neighbouring site, to the right or to the left with equal probability, with a rate dk only 
depending on the number of particles iVj = k, present on site %. As above, we consider 
a hopping rate such that dk ~ 1 + b/k at large k. The numerical simulations described 
below were performed with the choice of rate given by (14). 

4-1. Condensed phase (p > p c ) 

We begin with the determination of the dynamical exponent z. In regime II (k and t 
simultaneously large), let us assume the scaling form, inspired by eq. (21), 

/*(*) « W' 9 (^) • 
Then, at large time, the mean squared occupation 

oo 

^ = (*?(*)) = E*V*(0 (si) 

fc=l 

should behave as 

roc 

fi 2 ~t 1/z / duu 2 g{u). (32) 
Jo 

This quantity is easy to measure, and remains accurate even for relatively large times. 
More generally the moment of order n, p n = (N?(t)), scales as t^AXn-i). The behaviour 
(32) is very well observed in numerical simulations, with z = 3 (see figure 4 below). 
Figure 3 depicts the normalized scaling function ug(u)/(p — p c ) for times equal to 10 3 , 
3 10 3 , and 10 4 , with u = kt~ 1/3 . The data are for b = 4, and p = 20 p c = 10. 

We checked that the observed value of the dynamical exponent z = 3 was robust 
to changes in b. On the other hand, while for the mean-field case the scaling function 
g(u)/(p — p c ) was independent of the density p, in the one-dimensional case numerical 
simulations show dependence in the density. However this dependence can be absorbed 
in a rescaling of time: 

* — (P ~ Pc)t, (33) 

as demonstrated by figure 4. It depicts a plot of /^/(p — Pc) against ((p — p c )t) 1 / 3 , 
for p = 5, 10, 20. In the linear part of the curves, the slopes are all approximatively 
equal to 3, in very good agreement with the measured numerical values of the integral 
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0.2 




Figure 3. One-dimensional system: normalised scaling function ug{u)/(p — p c ) 
in the condensed phase for times equal to 10 3 , 310 3 , and 10 4 . (6 = 4, 
p = 20 Pc = 10.) 

in the right side of (32), divided by (p — p c ) 4 ^ 3 . Coming back to figure 3, plotting 
ug(u)(p — p c )~ 2 / 3 against k((p — pjt)" 1 / 3 for several densities leads to data collapse, 
confirming the rescaling of time by the density. 

We then determined the characteristics of regime I (k fixed, t large), for 6 = 4 and 
p = 10. We first observed that eqs.(18), (19) and (20) still hold. This was checked in 
a number of ways. For instance one can plot (/o//o,eq — 1) _1 and (d t — l)~ 1 against /i2, 
yielding straight lines of slopes respectively equal to v o I — —136 and A I ~ 59. This 
leads to A ~ 0.97, and vq — —2.2. These values are consistent with the relation (23), 
since the measured value of the integral / dugiu) ~ 1.7. The amplitude A also depends 
on p, with values equal respectively to 1.25, 0.97, and 0.76 (for p = 5, 10,20), which, 
when multiplied by (p — p c ) 1//3 , consistently with (33), collapse to a value ~ 2. 

4-2. Criticality (p = p c ) 

The critical case is harder to analyse numerically, because of the presence of large 
fluctuations in the system, and because the asymptotic regime is longer to attain. 
Guided by the mean-field analysis, we expect the dynamics to exhibit two scales of 
time, corresponding to regime I and II respectively. In the latter, let us assume the 
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Figure 4. One-dimensional system: rescaled mean squared occupation /^/(p — 
p c ) in the condensed phase, for three values of the density: p = 5, 10, 20, from 
bottom to top. (b = 4.) 



scaling behaviour, inspired from (24), 

fk(t) ~ fk, eq 9c(u) u = kr 1,z % 

where z c is the critical dynamical exponent. In order to determine the value of this 
exponent, we measured 

/•, I / 1 roo 

(34) 



Vb,; 



gMjfcn^fdWM) rdug c (u), 

k=0 Jk,eq J® 



for n = and n = 1, and where the notation reminds the fact that, at large values 
of k, fk,eq ~ k~ b . Figure 5 depicts a plot of v bQ for b = 3,4,5 against t 1 / 5 , which 
points towards z c = 5. Actually the best regression fit is obtained for z~ x ~ 0.215. We 
checked that u b>1 ~ (^,o) 2 - The ratio of these two quantities plotted against time gives 
an indication of how the asymptotic regime is attained. Figure 6 depicts a plot of the 
critical scaling function g c (u), against the scaling variable u = kt' 1 ^ 5 , for three different 
times (800, 2000, and 4000). 

In regime I the sum rules (1) and (2) yield, along the same lines as in section 3.2 
(see eqs. (25) and (26)), 
b-2 



(35) 



which generalizes (27). 




1 1 1 1 1 1 1 1 1 1 1 

1 2 3 4 5 6 



Figure 5. One-dimensional system: plot of Vbfl against t 1 ' 5 at criticality, for 
b = 3, 4, 5, with regression lines, (p = p c = 1/(6 — 2).) 

4-3. Universality 

As for the mean-field case, we compared the numerical results obtained with the choice 
of rates (14), (29), and (30). For instance, in the condensed phase, consider the mean 
squared occupation /i 2 , defined by equation (31). The global appearance of the set of 
three curves thus obtained is very similar to the corresponding set of curves obtained 
in the mean-field case, for the same choice of rates. For the latter, one observes weak 
discrepancies between the curves at shorter times, which however disappear at very large 
times, say 10 5 , obtained by numerical integration of the master equation (17). 

The same analysis can be done at criticality, for instance by considering the 
quantities v^q and \ defined in (34), with the same conclusions. In particular the 
ratio f6,i/(^6,o) 2 converges at large times to a universal amplitude ratio only dependent 
on b. 

We are thus lead to infer that, for the one-dimensional system, the scaling functions 
describing the asymptotic behaviour of the occupation probabilities fa for the 

mean-field system, universal quantities, independent of the detailed form of the hopping 
rate dk, provided that at large values of k, dk ~ 1 + bjk. However, in contrast with 
the mean-field case, in the condensed phase there exists a dependence on the density 
through a rescaling of time. 
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Figure 6. One-dimensional system: plot of the critical scaling function g c (u) 
for three different times (800, 2000, and 4000). (b = 4, p = p c = 0.5.) 



5. Discussion 

5.1. Universality: a summary, and the role of a bias 

Much emphasis was given in the present work to universal aspects of dynamics. 
Numerical simulations demonstrate the existence of scaling in the one- dimensional 
system, very akin to the scaling behaviour present at the mean-field level. The scaling 
regime is characterized by scaling functions, g and g c , in the condensed phase and at 
criticality respectively, which do not depend on the detailed form of the rate dk, provided 
that at large values of the occupation variable, dk ~ 1 + b/k. They only depend on the 
parameter b, and, in the condensed phase, on the density, in contrast with the mean- 
field case. However a rescaling of time by the density, and correspondingly of g, leads 
to excellent data collapse. This universality extends to the case of the mean-field zeta 
urn model, though the definition of the model is different from that of the zero-range 
processes studied here. The dynamical exponents z and z c are universal in a stronger 
sense since they neither depend on b, nor on the density. They only depend on the 
dimensionality of space, with the mean-field geometry corresponding to the limiting 
case d — oo. 

The observed differences between the values of the dynamical exponents z and z c , in 
the mean-field geometry on one hand, and in the one-dimensional case on the other hand, 
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give a measure of the role of the diffusion of particles from sites to sites. Indeed, besides 
the fact, independent of the geometry of the system, that the occupation numbers iVj 
are fluctuating variables by the very definition of the process, for the one-dimensional 
system a new phenomenon appears, which is the spatial diffusion of particles along the 
line. This induces correlations between neighbouring boxes (sites). 




Figure 7. One-dimensional system: asymmetric ZRP (p = 1). 

Left: rescaled mean squared occupation /U2/G ~~ Pc) against i 1 / 2 in the 

condensed phase, for p = 5, 10, 20, from bottom to top. (b = 4.) 

Right: 1^,0 against t 1 / 3 at criticality, for b = 3, 4, with regression lines. 

(p = Pc = 1/(6 - 2).) 

To gain more understanding of this point, we introduced a bias in the move of a 
particle out of a box. A particle now hops to the right with probability p and to the left 
with probability q — 1 — p. This does not change the stationary state. For any positive 
value of the bias p— q we found z = 2, and z c ~ 3 (see figure 7). Therefore the symmetric 
and asymmetric one- dimensional ZRP belong to two different universality classes. In 
the extreme case where p = 1, corresponding to the totally asymmetric ZRP, there is no 
diffusion left in the motion of a particle along the line. One could therefore expect that 
the totally asymmetric ZRP should resemble, at least to some extent, to the mean-field 
ZRP. It indeed turns out that, in the condensed phase, the scaling function g{u) of the 
former is quantitatively very close to the scaling function of the latter. Moreover, for 
p — 1, the universality of the scaling function g{u) with respect to density is restored, 
without rescaling of time. (This point was also checked for the value p = 3/4.) However 
the fact that, seemingly, z c ~ 3 is a caveat that the influence of a bias is more subtle. 
As a first step, it would be interesting to have a better understanding of the role of the 
bias on the decay of correlations. 

A summary of the different cases studied in this work is presented in table 1. 
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Table 1. Values of the dynamical exponents for three classes of universality: mean- 
field geometry (d = oo), one-dimensional asymmetric (d = l,p / |), onc-dimensional 
symmetric (d = l,p = j). 





d = oo d=l,p^\ d = 




z 

Zc 


2 2 3 





5.2. Future work 

Finding an analytical argument in favour of the values of the exponents appearing in 
the table is the most immediate challenge. Other immediate points of interest are the 
determination of the upper critical dimension for the zero-range processes under study, 
or the explanation of the role of the density in the one-dimensional condensed phase. 

The value z = 3 is reminiscent of the value of the dynamical exponent for the 
growth of domains in the dynamics of an Ising chain under Kawasaki dynamics, in the 
scaling regime of low temperature [25]. It is interesting to note that, for the Kawasaki 
chain, adding a bias induces a change in the dynamic exponent from z = 3 to z = 2 
[26, 27]. Moreover it is simple to see that a mean-field version of the model leads to 
z = 2. Think to the low-temperature regime, where the detachment of a spin from a 
domain is a rare event. Then whenever such an event occurs, one of the spins which 
are released (e.g. the (+) spin) is immediately transferred to a random location inside 
a domain of the opposite sign. By diffusion it will shortly stick to the nearest domain 
of its sign, thus increasing the size of this domain by one unit, to the expense of the 
initial domain, far away in the system. It is clear from this analysis that the mechanism 
leading to the value z = 3 in the low-temperature coarsening regime of the Kawasaki 
chain [25] no longer works, hence z — 2. Other mean-field versions in the same spirit 
can be invented, with the same result regarding the value of z. 

Trying to draw a closer analogy, beyond these simple observations, between the 
dynamics of the Ising chain under Kawasaki dynamics and the dynamics of the one- 
dimensional ZRP deserves further analysis. 

Another point of important interest concerns non-stationary aspects of the 
dynamics. Non stationarity is demonstrated by the behaviour of two-time quantities 
such as two-time autocorrelation, response and fluctuation-dissipation ratio. At the 
mean-field level, guided once more by the analysis of the dynamics of the zeta urn 
model [22], we can show that the formalism developed for this model holds for the ZRP 
studied in the present work. This, as well as the study of the behaviour of two-time 
quantities for the one-dimensional ZRP, will be the subject of another work [24]. 

The last challenging point concerns the possible relevance of the present work for the 
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description of the class of driven particle systems exhibiting a condensation transition, 
such as those studied in [19]. The question posed is whether these particle systems can 
be described during their temporal evolution, and not only at stationarity, in terms of 
the zero-range processes studied in the present work [20]. As a first step one can use 
the present study as a template, and ask the following questions [19]. In the condensed 
phase, or at criticality, is there a scaling regime, and if so, what are the values of the 
dynamical exponents? What is the role of the density? What is the influence of a bias? 
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